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Abstract 

I discuss the theoretical uncertainties in the extraction of | V c b\ 
from a measurement of the B — > D*£ v decay rate close to zero recoil. 
In particular, I combine previous estimates of the l/mg corrections to 
the normalization of the hadronic form factor at zero recoil with sum 
rules derived by Shifman et al. to obtain a new prediction with less 
uncertainty. I also give a prediction for the slope of the form factor 
£(w) at zero recoil: g 2 = 0.7±0.2. Using the most recent experimental 
results, I obtain the model-independent value | V c b\ = 0.0395 ± 0.0030. 
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1 Introduction 



With the discovery of heavy quark symmetry (for a review see Ref. [|l| and 
references therein), it has become clear that the study of exclusive semilep- 
tonic B — > D*£ v decays close to zero recoil allows for a reliable determination 
of the Cabibbo-Kobayashi-Maskawa matrix element V c b, which is free, to a 
large extent, of hadronic uncertainties [0-0]. Model dependence enters this 
analysis only at the level of power corrections, which are suppressed by a 
factor of at least (AQ CD /m c ) 2 . These corrections can be investigated in a 
systematic way using the heavy quark effective theory ||. They are found 
to be small, of the order of a few per cent. 

Until recently, this method to determine | V c b\ was limited by large ex- 
perimental uncertainties of about 15-20%, which were much larger than the 
theoretical uncertainties in the analysis of symmetry-breaking corrections. 
However, three collaborations have now presented results of higher precision 
It is thus important to reconsider the status of the theoretical analy- 
sis, even more so since the original calculation of power corrections in Ref. 
has become the subject of controversy 
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The differential decay rate for the process B — > D*£ v is given by [1 j 



dT 

dw 



487T 3 

1 + 
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v cb \ 2 v 2 A i 2 H 



where 



m\ + mf>. - q 2 , „ 

w = v B ■ v D * = (2) 

zrriBmD* 

denotes the product of the meson velocities. I have factorized the hadronic 
form factor for this decay into a short-distance coefficient tja and a func- 
tion £(u>), which contains the long-distance hadronic dynamics. Apart from 
corrections of order 1/rriQ, this function coincides with the Isgur-Wise form 
factor H pj| . Luke's theorem determines the normalization of £,{w) at zero 
recoil (w = 1) up to second-order power corrections 0, 12]: 



£(i) 



/m 2 



(3) 
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The strategy is to obtain the product | Va,\r]A^(w) from a measurement of 
the differential decay rate, and to extrapolate it to w — 1 to extract 

I V cb \ r) A (1 + 6 1/m ») = | V cb \ {l + 0(a s , l/m 2 Q )) . (4) 

The task of theorists is to provide a reliable calculation of rjA and Si/ m 2 in 
order to turn this measurement into a precise determination of | V cb \. 



2 Calculation of t]a 

The short-distance coefficient T)a takes into account a finite renormalization 
of the axial vector current in the region m b > fi > m c . Its calculation is 
a straightforward application of QCD perturbation theory. At the one-loop 
order, one finds [fj, [13| 

a s fm b + m c m b 8\ 

VA = 1 + — In • 5 

7r \m b — m c m c 3/ 

The scale of the running coupling constant is not determined at this order. 
Choosing a s between a s {m b ) ~ 0.20 and a s {m c ) ~ 0.32, and using m c /m b = 
0.30 ± 0.05, one obtains values in the range 0.95 < t/a < 0.98. The scale 
ambiguity leads to an uncertainty of order At]a ~ [(a s /ir) \n(m b /m c )} 2 ~ 2%. 

The calculation can be improved by using the renormalization group to 
resum the leading and next-to-leading logarithms of the type [a s ln(m b /m c )} n , 
a s [a s \n(m b /m c )] n , and {m c /m b )[a s lia(m b /m c )] n to all orders in perturbation 
theory [14j]-[T7]]. A consistent scheme for a next-to-leading-order calculation 
of rjA has been developed in Ref . [H| . The result is 



rj A = x 6/25 |l + 1.561 



a s {m c ) - a s {m b ) _ 8a s (m c ) 

it 3n 
, 25 _ 14 x _ 9/25 1_ x _ 12/25 _8_ ■ 

m b V54 27 18 25 

| 2a s (m) m 2 c \ n mb \ ( 6 ) 
7r m b {m b — m c ) m c j ' 

where x = a s (m c )/a s (mb), and > m > m c . The numerical result is very 
stable under changes of the input parameters. For Aj^g = (0.25 ± 0.05) GeV 
(for 4 flavours) and m c /m b = 0.30 ± 0.05, one obtains t/a = 0.985 ± 0.006. 
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The uncertainty arising from next-to-next-to-leading corrections is of order 
At]a ~ (a s /ir) 2 ~ 1%. Taking this into account, I think it is conservative to 
increase the error by a factor 2.5 and quote 

T]a = 0.985 i 0.015 . (7) 

3 Anatomy of #i/ m 2 

Hadronic uncertainties enter the determination of | V c b\ at the level of second- 
order power corrections, which are expected to be of order (Aqcd/^c) 2 ~ 3%. 
For a precision measurement, it is important to understand the structure of 
these corrections in detail. Falk and myself have derived the exact expression 

§ 

which depends upon five hadronic parameters: tp and £y parametrize the 
deficit in the "wave-function overlap" between b- and c-flavoured pseudoscalar 
(P) and vector ( V) mesons, for instance 

(DW | c t HBW)Kl _(J___i.)% p . (9) 

The parameter — Ai = (pq) is proportional to the kinetic energy of the heavy 
quark inside a heavy meson, A2 = (m 2 B * — m 2 B )/A ~ 0.12 GeV 2 is determined 
by the vector-pseudoscalar mass splitting, and A^a parametrizes certain ma- 
trix elements containing double insertions of the chromo-magnetic operator. 
With the exception of A2, estimates of these parameters are model-dependent. 
In Ref. ||, we made the simplifying assumptions that £ P = £ v , and that the 
corrections represented by \ G i are negligible. The latter one is based on the 
observation that these corrections involve a double insertion of an operator 
that breaks the heavy quark spin symmetry. Using then reasonable values 
such as l P = i v = (0.35 ± 0.15) GeV 2 and -A x = (0.25 ± 0.20) GeV 2 , one 
obtains 6i/ m 2 = —(2.4 ± 1.3)%. Here and in the following, I take m b = 4.80 
GeV and m c = 1.45 GeV for the heavy quark masses. In Ref. ^ ne error 
in the estimate of 5i/ m a has been increased to ±4% in order to account for 
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the model dependence and higher-order corrections. A very similar result, 



—5% < 8x/ m 2 < 0, has been obtained by Mannel . 

Recently, Shifman et al. have suggested an alternative approach to obtain 
an estimate of 5i/ m 2, which is based on bounds derived using sum rules and 



the operator product expansion [10|. These bounds imply the inequalities 



^p>i(-Ai-3A 2 )=^ in >0 ! 
V>^(-Ai + A 2 )^ in >2A 2 



/ 1 1 \/£^ m £^ n \ 1 /4 . n . 

<W 2 < — — H Ai + 2A 2 

\2m c 2m^J \2m c zm^J Amcirib V3 

< -7^% - -2-9%. (10) 
2mi 



The upper bound for 5i/ m a implies that tja£(1) < 0.956. In Ref. fllPl , this 
number is quoted as 0.94. In the same reference, the authors give an "edu- 
cated guess" ^£(1) = 0.89 ± 0.03 corresponding to Sy m 2 = —(9.6 ± 3.0)%. 
However, the arguments presented to support this guess are not very rigorous. 
It is possible to combine the above approaches to reduce the theoretical 



uncertainty in the estimate of Su m 2 [p0[ . The idea is to use the sum rules 



to constrain the hadronic parameters in (|5p in a threefold way: (i) The first 
relation in (pTOj) implies that 

- Ai > 3A 2 ~ 0.36 GeV 2 , (11) 

excluding some of the values for the parameter Ai used in previous analyses. 

(ii) Comparing the third relation in (|T0| ) with (H) in the limit mj = m c , one 
finds that 

A G 2>0. (12) 

(iii) Finally, £p and £y are correlated in such a way that £y > tp if A G 2 is 
not too large. To illustrate this last point, let me define new parameters 



i=\(lv + lp), D= l -(£ v -£ P )-\ 2 . (13) 



In terms of these, 



6i/m2 = - ~ 1 ~ (ii ~ if J (A2 + D) 
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Figure 1: Allowed regions for Su m 2 as a function of X G 2 for the two cases 
1= 0.2 GeV 2 (solid) and 0.4 GeV 2 (dashed). 



+ 
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4m c mf, \3 



Ai + 2A2 — Xg 2 ■ 



(14) 



Using the inequalities fllP]) , one can show that £ > |(— Ai — A 2 ) and — -D max < 

d < J o max , with m 



s 



< S < \ G 2/2 



A G2 5-A 2 2 /4; S> A G2 /2. 



(15) 



where S = I + (Ai + A 2 )/2. There are thus three effects, which decrease 
5i/ m 2 with respect to the estimate given in Ref. 0: a large value of (— Ai), a 
positive value of A G 2, and the fact that for small A G 2 the difference (£y — £p) 
is centred around 2A 2 (i.e. D is centred around 0). 

In evaluating (|T4]), I will take — Ai = 0.4 GeV 2 , which is consistent with 
the bound in (pTTj) and with the value — Ai = (0.5 ± 0.1) GeV 2 obtained from 
QCD sum rules [2l]|. Varying — Ai in the range between 0.36 and 0.5 GeV 2 
does not alter the results very much. The main uncertainty comes from 
the unknown values of the parameters I and A G 2. As a guideline, one may 
employ the constituent quark model of Isgur et al. |22|], in which one uses 
non-relativistic harmonic oscillator wave functions for the ground-state heavy 



mesons, for instance ipB{ r ) 



exp(-|/iu;r 2 ) 



where /1 = (l/m q + l/m&) 
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is the reduced mass. One then obtains I = |m 2 ~ 0.2 GeV 2 , where I use 
m q ~ 0.5 GeV for the light constituent quark mass. However, this estimate of 
I is probably somewhat too low. Lattice studies of heavy-light wave functions 
suggest an exponential behaviour of the form ^(r) ~ exp(— k\it) [p3| |, which 
leads to I = |m 2 ~ 0.4 GeV 2 . Values much larger than this are unlikely, since 
I use a rather large constituent quark mass m q . In fact, adopting the point of 
view that the sum rules for £ P and ty are saturated to approximately 50% by 



the ground-state contribution |10| , one would expect i ~ (— Ai — A2) — 0.28 
GeV 2 , which seems a very reasonable value to me. In Fig. [l|, I show the 
allowed regions for 5i/ m a as a function of Aq2 for two values of I. I think it 
is reasonable to assume that \qi is of a magnitude similar to A2 or smaller. 
Thus, I conclude that for all reasonable choices of parameters the results are 
in the range 

- 8% < 6 1/m 2 < -3% , (16) 

which is consistent with the previous estimates in Refs. ||, |K| at the la 
level. A more precise determination of the parameter I would help to reduce 
the uncertainty in this number. 



4 Prediction for the slope parameter g 2 

In the extrapolation of the differential decay rate (JI|) to zero recoil, the slope 
of the function £(u>) close to w = 1 plays an important role. One defines a 
parameter g 2 by 

£H=£(l){l-£ 2 ( W -l) + ...}. (17) 

It is important to distinguish q 2 from the slope parameter q 2 of the Isgur- 
Wise function. They differ by corrections that break the heavy quark sym- 
metry. Whereas the slope of the Isgur-Wise function is a universal, mass- 
independent parameter, the slope of the physical form factor depends on 
logarithms and inverse powers of the heavy quark masses. On the other 
hand, g 2 is an observable quantity, while the value of g 2 depends on the 
renormalization scheme. To illustrate this last point, let me neglect for the 
moment 1/rriQ corrections and work in the leading logarithmic approxima- 
tion. Then the relation between the physical slope parameter q 2 and the 
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slope parameter g 2 (fi) of the regularized Isgur-Wise function is 



^9 9/ n 16 , Oi s {m) 9 16 , , < 

e = ^-ir ln ^ se '-si "»«•(•»)■ < 18 > 

where /i is the renormalization scale, and m is an undetermined (at this 
order) scale between m;, and ra c . The last equation can be used to define the 
/^-independent slope of the renormalized Isgur-Wise function (see Ref. [[18] 
for the generalization of this definition to next-to- leading order). If next-to- 
leading logarithmic corrections are taken into account, the scale ambiguity 
related to the choice of m is resolved, and one obtains 

£ 2 = £ 2 + (0.14 ± 0.02) + 0(l/m Q ). (19) 

The 1 / itiq corrections to this relation have been investigated and are found 
to be negative. However, any such estimate is model-dependent and thus has 
a large theoretical uncertainty. The result is q 2 ~ q 2 ± 0.2 |2(J. Predictions 
for the renormalized slope parameter q 2 are available from QCD sum rules 
including a next-to- leading-order renormalization-group improvement. One 
obtains q 2 ~ 0.7 ± 0.1 H |§ H). I thus expect 



f = 0.7 ±0.2. (20) 



5 Summary 



Using the updated values t/a = 0.985 ± 0.015 and <5i/ m 2 = —(5.5 ± 2.5)%, I 
obtain for the normalization of the hadronic form factor at zero recoil: 



77^(1) = 0.93 ±0.03 



(21) 



Three experiments have recently presented new measurements of the product 
I K*| 7?a£(1)- When rescaled using the new lifetime values t b o = (1.61 ± 
0.08) ps and t b + = (1.65 ± 0.07) ps [26], the results obtained from a linear 
fit to the data are 



0.0347 ±0.0019 ±0.0020 
V cb \ tia£Q) = { 0.0364 ± 0.0042 ± 0.0031 
0.0385 ±0.0043 ±0.0028 



Ref. [6], 
Ref- 0, 
Ref. B, 



(22) 
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where the first error is statistical and the second systematic. I will follow 



the suggestion of Ref. [g7j an d add 0.001 ± 0.001 to these values to account 
for the curvature of the function Taking the weighted average of the 

experimental results and using the theoretical prediction (|21|), I then obtain 

| V cb \ = 0.0395 ± 0.0027 (exp) ± 0.0013 (th) = 0.0395 ± 0.0030 , (23) 

which corresponds to a model- independent measurement of | V c t,\ with 7% 
accuracy. This is by far the most accurate determination to date. 

Neglecting 1/mq corrections, I have related the physical slope parameter 
g 2 to the slope of the Isgur-Wise function and obtain the prediction g 2 = 
0.7 ± 0.2. It compares well with the average value observed by experiments, 
which is g 2 = 0.87 ± 0.12 §-[g. 



References 

[1] M. Neubert, SLAC preprint SLAC-PUB-6263 (1993), to appear in Phys. 
Rep. 

[2] M.B. Voloshin and M.A. Shifman, Yad. Fiz. 47, 801 (1988) [Sov. J. Nucl. 
Phys. 47, 511 (1988)]. 

[3] N. Isgur and M.B. Wise, Phys. Lett. B 232, 113 (1989); 237, 527 (1990). 

[4] M. Neubert, Phys. Lett. B 264, 455 (1991). 

[5] H. Georgi, Phys. Lett. B 240, 447 (1990). 

[6] T. Browder (CLEO Collaboration), talk presented at the 27th Interna- 
tional Conference on High Energy Physics (ICHEP94), Glasgow, Scot- 
land, July 1994; B. Barish et al. (CLEO Collaboration), Cornell preprint 
CLNS 94/1285, submitted to ICHEP94. 

[7] I. Scott (ALEPH Collaboration), talk presented at ICHEP94; D. De- 
camp et al. (ALEPH Collaboration), paper gls0605, submitted to 
ICHEP94. 

[8] H. Albrecht et al. (ARGUS Collaboration), Z. Phys. C 57, 533 (1993). 



8 



[9] A.F. Falk and M. Neubert, Phys. Rev. D 47, 2965 and 2982 (1993). 

[10] M. Shifman, N.G. Uraltsev, and A. Vainshtein, Minnesota preprint TPI- 
MINN-94/13-T (1994); I. Bigi, M. Shifman, N.G. Uraltsev, and A. Vain- 
shtein, Minnesota preprint TPI-MINN-94/12-T (1994). 

[11] A.F. Falk, H. Georgi, B. Grinstein, and M.B. Wise, Nucl. Phys. B 343, 
1 (1990). 

[12] M.E. Luke, Phys. Lett. B 252, 447 (1990). 
[13] M. Neubert, Nucl. Phys. B 371, 149 (1992). 

[14] H.D. Politzer and M.B. Wise, Phys. Lett. B 206, 681 (1988); 208, 504 

(1988) . 

[15] X. Ji and M.J. Musolf, Phys. Lett. B 257, 409 (1991). 

[16] D.J. Broadhurst and A.G. Grozin, Phys. Lett. B 267, 105 (1991). 

[17] A.F. Falk and B. Grinstein, Phys. Lett. B 247, 406 (1990). 

[18] M. Neubert, Phys. Rev. D 46, 2212 (1992). 

[19] T. Mannel, Phys. Rev. D 50, 428 (1994). 

[20] M. Neubert, CERN preprint CERN-TH. 7395/94 (1994), to appear in 
Phys. Lett. B. 

[21] P. Ball and V.M. Braun, Phys. Rev. D 49, 2472 (1994), and private 
communication. 

[22] N. Isgur, D. Scora, B. Grinstein, and M.B. Wise, Phys. Rev. D 39, 799 

(1989) . 

[23] A. Duncan, E. Eichten, and H. Thacker, Phys. Lett. B 303, 109 (1993). 
[24] E. Bagan, P. Ball, and P. Gosdzinsky, Phys. Lett. B 301, 249 (1993). 
[25] M. Neubert, Phys. Rev. D 47, 4063 (1993). 
[26] F. DeJongh, talk presented at ICHEP94. 
[27] R. Patterson, talk presented at ICHEP94. 



9 



